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Abstract 

As previously observed the massless limit of massive gravity leads 
to a modification of general relativity. Here we study spherically sym- 
metric solutions of the modified field equations which contain normal 
matter together with a dark energy density. If the dark density profile 
is assumed to be known, the whole problem is reduced to a linear first 
order differential equation which can be solved by quadratures. 
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1 Introduction 



In a previous paper with the same title [1] we have made the thousand 
and first proposal to explain dark matter density. In contrast to the other 
thousand proposals we came to this subject by accident, that means not by 
looking for an explanation of dark matter. As the title indicates we have 
studied massive quantum gravity which is gravity with a massive graviton. 
To have a massive spin-2 gauge theory the so-called vector graviton field v\ 
is indispensable. The crucial observation was that in the limit m — > of 
vanishing graviton mass this field v\ does not decouple from the symmetric 
tensor field h^ u which is equivalent to Einstein's g^ v . In the classical limit v\ 
acts as four scalar fields v n . As a consequence the massless limit of massive 
gravity is different from general relativity. There remains the additional 
coupling to the four (now massless) scalar fields v n . In the resulting modified 
Einstein's equations this gives additional terms in the energy - momentum 
tensor. In the 00-component the new contribution looks as if it comes from a 
dark matter density, but there are also peculiar modifications of the pressure 
components in the jj-equations. 

The paper is organized as follows. In the next section we review the 
modified general relativistic equations, for their derivation we refer to [1]. 
We then consider static spherically symmetric solutions including normal 
matter which is described by a mass density q(r) and an isotropic pressure 
p(r). By expanding the solution for large distance r we get already an 
important result: the dark density profile must decrease as 1/r 4 as already 
found in [1] for the solution without matter. This tail contradicts the widely 
discussed Navarro - Frenk - White (NFW) profile [2] [3] 

P{r) = (r/r a )(l + (r/r a ))2 

which has a 1/r 3 tail. However, this profile must be modified for large r 
anyway to get a finite total mass for the dark halo. A second conclusion can 
be drawn from the expansion around r = 0. In general relativity there exists 
the inner Schwarzschild solution which is finite in all quantities including 
mass density and pressure. This is of course also a solution of our modified 
theory with vanishing dark density. We can now test whether there is a 
corresponding finite solution with a dark density different from zero. The 
answer is no. This may explain the fact that a single star like the sun which 
is described by this finite solution does not have a dark halo. In order to 
describe dark halos we have to study solutions with some singularity for 
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small r. In [1] we have investigated the vaccum solutions, here we start the 
analysis of solutions with normal matter. 

The modified general relativity has the following nice property. To spec- 
ify a solution one quantity must be given which usually is the mass density 
q(r) of normal matter. In this case one has to solve non-linear differential 
equations. Alternatively, if the dark profile p(r) is considered to be known, 
the remaining differential equation is linear and first order. It can be solved 
by quadratures. This is a machine which produces solutions at low cost. In 
this way we construct a simple singular solution in section 3. This solution 
with a dark halo has a singularity at some finite radius r s . There is no 
horizon at r > r s so that the singularity is a naked one. 



2 Massless limit of massive gravity 

The basic classical field equations which follow from massive gravity in the 
limit of vanishing graviton mass [1] are the modified Einstein equations 

R^u-^g^R-^^-t^ = ^^(-d^v n d u v n ++^g^g a Pd a v n dpv n y (2.1) 

together with the Laplace-Beltrami wave equation for the four scalar fields 
v n ,n = 0,1,2,3 

d a (^gg at5 d p v n ) = 0. (2.2) 

As discussed in [1] the Latin index n of v n is raised and lowered with the 
Minkowski tensor diag(l, —1,-1,-1), in contrast to the Greek indices which 
are changed with the metric tensor g^ v . Since upon v n we act by partial 
derivatives not by covariant ones, this means that the v n are four scalar 
fields, indeed. The new contribution on the right-hand side of (2.1) is the 
possible origin of the dark density. 

We want to study static spherically symmetric solutions of these field 
equations. As in [1] we write the metric as 

ds 2 = e u c 2 dt 2 - e x dr 2 - r 2 ((tt) 2 + sin 2 tidip 2 ) (2.3) 

where v and A are functions of r only. We take the coordinates x° = ct, 
x 1 = r, x 2 = i?, x 3 = if such that 

5oo = e", gu = -e x 

922 = ~r 2 , 533 = ~r 2 sin 2 ■d (2.4) 
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and zero otherwise. The components with upper indices are the inverse of 
this. The determinant comes out to be 

g = detg^ = -e u+x r 4 sin 2 (2.5) 

The energy momentum tensor in (2.1) is assumed in the simple form 

t£ = diag^c 2 , -p, -p, -p) (2.6) 

where q(r) is the ordinary mass density and p(r) is an isotropic pressure. 
Now the following modified radial Einstein equations must be solved 



« (7 - I2J+I2 = -TfoWc + Mr)) (2-7) 



\(v' 1 \ 1 8ttG ,p(r) ... . 

e - + 12 -12 = -^-(^r " w oW) (2.8) 



r r 2 J r 2 c° c 



x( u" u' 2 v'\' X' v'\ 8nG /P (r) , 

e (Y + T " — " 2; + 2F)= ~ ^ + Wo(r)) ' (2 - 9) 

where t^o is the contribution of the four scalar fields v n . As discussed in [1] 
in the spherically symmetric case this is of the form 

First by suitable combination we simplify the equations. Adding (2.7) 
to (2.8) we get 

e —(\l + V >)=g( qC+ l) (2.11) 

r c 

where 

g =—. (2.12) 
Eliminating p from (2.8) and (2.9) we obtain 

z/' = 2e A (-- + 2^ +— - — + + -. 2.13 

\ r z c J 2 2 r r A 

Next we differentiate (2.8) with respect to r and use (2.13): 

_S4) = _.-»£ ( V + l0 + (2.14) 
c / 2r cr 
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Substituting (2.11) inhere we finally arrive at 

p' v' ( p 2 



+ (2.15) 



C 2 2 v c 

where (2.10) has been used. This differential equation for the pressure will 
be used instead of the second order equation (2.9) in the following. 

To study the solution for large r where it should approach flat space, we 
set up an expansion in powers of 1/r: 

e- A = i-7 + f§ + - <"<0 

„ = £ + £ + (2.17) 
,-^ + ^ + ... (2.18) 
!>=jS + jj + - (2-19) 



We find 



oi = -6i (2.20) 

a2 = 5(94 + Pi), Pi = — (2-21) 

c 

94 = -(4& 2 + fo?)-2pi (2.22) 



2 = -77^(94 + 2pi) = -^(462 + 6f). (2.23) 
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To (94 + 2pl) = -io 5 

If the dark density u>o(j") is given (2.10), i.e. po and i/(r) are known, we 
have a unique solution. The equation of state of the normal matter is then 
determined. 

In ordinary general relativity there exists a finite inner solution which 
can be expanded around r = in the form 

e" A = 1 + air + a 2 r 2 + a 3 r 3 + . . . 

v = b + b 1 + b 2 r 2 + b 3 r 3 + ... (2.24) 

qc = q + qir + q 2 r 2 + ... 

P 2 

- = Po +Pir + p 2 r + 

c 
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Substituting this into (2.7-9) with u>o = we find 

g g 

ai = 0, a 2 = -gtfo, «3 = -^Qi (2.25) 

bi=0, 62 = |(P0 + f), ^3 = ^1 (2.26) 

Pi=0, P2 = -|(po + f )(po + ?o) (2-27) 

P3 = -gPO + Y^909l)- 

Now we can check whether this finite solution has a counterpart with non- 
vanishing wq. Expanding (2.10) 



.4 o ' • • • 



we easily see from (2.7) that po must vanish. A solution finite at r = 
cannot have a dark halo. Therefore we have to look for singular solutions. 

3 A singular solution 

To construct a singular solution with non- vanishing dark density Wo(r) we 
proceed as follows. We eliminate p(r) in the two equations (2.8) (2.9): 



" Ae (T + 27) =e (-Y-T + 27 + ^)- 

-^ + 2gw . (3.1) 



Introducing 

V = e~ x (3.2) 
we have obtained a linear first order differential equation for y(r): 

h(r)y'{r) = f 2 (r)y + h(r) (3.3) 

where 

i/" z/ 2 i/' 1 
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f 3 = -^ + 2gw . (3.6) 



Considering the dark density wo(r) (2.10) as given we can compute 

v{r) = log po — log wo — 4 log r. (3-7) 
This enables us to express the coefficients /i, /b, in terms of w$: 



i< 3w 2 5 w 
2wo 4«)q 2 rwg ^ 2 



A ^ T^r ^ 3 (3.9) 



/3 = -^2+2^o- (3.10) 

The remaining equation (2.7) then determines the mass density q(r) of the 
normal matter. 

The linear equation (3.3) can be simply solved by quadratures, so we 
have full control of the solution. The general solution y(r) is the sum of a 
particular solution y\{r) of the inhomogeneous equation (3.3) plus a solution 
yo(r) of the homogeneous equation. We will soon realize that the solution 
is uniquely fixed by requiring that it approaches flat space for r — > oo. To 
see this we set up an expansion of the form 

» = + £ + - P.") 

In agreement with the expansion (2.16-19) above the dark density (2.10) 
must start as follows 

^0 = ^4 +^ + ^ + --- ( 3 - 12 ) 

A 1/r 3 tail is in contradiction with the flat space asymptotic y — > 1 for 
r -► 00. From (3.12) and (3.8) we find 

2r 4/9or 2 V2p 4pg/ r 6 

and similarly 

/2= ^ + ^ +( i*_?4)^ + ... (3 , 4 ) 
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/3 = -^ + 2 9 (f2 + * + * + ...). (3.15) 

Substituting all this into (3.3) we can determine the coefficients in (3.11). 
The result is 

d x ,2<% 2d 2 \ 1 , . 

y = i + (-9 5PoH + --- ( 3 - 16 ) 



Po?" v Po Po 1 r z 

To carry out all integrations in terms of elementary functions we choose 
the following simple dark density profile 

wo = — - — -7. (3.17) 

l + (r/r s ) 4 v ' 

It contains two parameters w\ and r s as most phenomenological dark matter 
profiles in the literature and has the correct asymptotic behavior. First we 
want to calculate the solution yo(r) of the homogeneous equation 

hVo = f2W, (3-18) 

which is given by 

r 

yo{r) = exp J j j-dx. (3.19) 
Now from (3.8) and (3.17) we get 

4 4 

h = ^fr y (3-20) 

In view of the integral in (3.19) it is convenient to use the original definitions 
(3.4) (3.5) which yield 

fi fi r r 2 f ± 

In our special case (3.17) this is equal to 

rrA-'-Ay < 3 - 22 > 

This can be easily integrated and exponentiated using (2.10) 

— v 4^0 

e = r — . 

Po 
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In this way we find the following solution (3.19) 



14 f 4 + r i\ 

y (r) = const. (r L,4)4 • ( 3 - 23 ) 



The solution of the linear homogeneous equation contains a free prefactor, 
of course. For r — ► oo yo( r ) increases as r 2 which is in conflict with the 
asymptotic flatness y ~ const. Consequently, the prefactor of yo has to be 
zero so that we get a unique solution of our problem as in the expansion 
(2.16) in the last section. 

From yo(r) on finds the solution yi(r) of the inhomogeneous equation by 
the well-known method of variation of the constants. One makes the ansatz 

j/i (r) =a(r)y (r), (3.24) 

where yo is given by (3.23) with prefactor 1. Substituting this into (3.3) one 
obtains a(r) by another integration: 

air) = J /3(g) dx, (3.25) 

where a constant of integration, i.e. the lower limit of the integral is free. 
The integration of the rational function (3.25) is elementary. The final result 
is 

- J_ _ r l 3 zl _ r s 12 

+4^4-2^(1 + -_-| + --!,]. (3.26) 

Substituting this back into (3.24) we get the desired unique solution which 
we denote by y(r) again: 

y( r ) = _ll±jj_f r 12 _ 4 8 + 3 8 4 _ I 12 

yy > (- r 4_ r 4^4\' 1 s> T 5's' 7 s 



+gwi 



-8r 14 log(l + T -f) + 7rfr 10 - 2r s 8 r 6 + ^r s 12 r 2 ] }. (3.27) 

The solution y{r) has the correct asymptotic behavior y — > 1 for r — > oo. 
However, there is a singularity at r = r s . We ask whether this singularity 
is hidden by a horizon. A horizon corresponds to y = exp(— A) = 0. So we 
look for a zero of the curly bracket in (3.27). If we insert realistic numbers 
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for the dark matter density w\ as discussed in the next section, we find that 
the second line in (3.27) is completely neglegible compared to the first one 
due to the smallness of g. The polynomial 

r 12_ V 3 8r 4_ 1 12 

5 7 

has only one real zero tq < r s . Consequently the singularity can be seen 
from the outside, that means it is naked. 

Finally the mass density q(r) of normal matter is obtained from (2.7) 
according to 

^) = - g (- V j-^ + ^)-w . (3.28) 
This also becomes singular at v — r s . 



4 Comparison with observations 

At first side it is questionable whether our solution with a spherically sym- 
metric density distribution q(r) of the normal matter can be compared to 
observations of spiral galaxies. However, since the dark halo is approxi- 
mately spherically symmetric, at least the gas component in the outer part 
of the halo can be assumed to be also spherically symmetric [4]. If we chose 
a galaxy which is dark matter dominated where the normal matter is only 
10 percent of the total mass, say, then a spherically symmetric description is 
not a bad approximation. The spiral galaxy M33 is a good example for our 
purpose. Edvige Corbelli [5] has obtained good fits to the data which give 
(via rotation curves) the density profiles of dark and visible matter. The 
data are compatible with very different forms of the dark density profile. 
Beside the NFW profile (1.1) Corbelli considers the isothermal profile 

*»W <4I) 

and the Burkert profile [6] 

Pb 

P{r) = (l + r/r B )(l + (r/r B ) 2 )' ^ 

which both have a central constant - density core like our profile (3.17). The 
data are not much constraining the profile in the central region, and in the 
tail region r > 17 kpc there are no data. Consequently it is well possible to 
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fit our profile (3.17) to the best fit of Corbelli between, say, 4 and 14 kpc. 
However, the resulting singular radius r s comes out to be around 9 kpc. 
This is too big. We could except a value r s < 0.5 kpc which is the radius of 
the bright "nucleus" [5] , where something unknown is going on. 

A similar conclusion can be drawn from the analysis of the following 
3-parameter profile 

w ° = ~YT\ , 2V ( 4 -3 

r z (l — a±r + d2r z ) 

Then the solution of the homogeneous equation (3.18) is given by 

, , r- 2 - 16a2 / a Ul-a 1 r + a 2 r 2 ) / 8\ 
Voir) = 72 —, 2 expl . (4.4) 

This produces a singularity at r s = a\/2a2- Again fitting (4.3) to Cor- 
belli's result for the dark matter profile we get r s around 10 kpc which is 
unacceptable. 

The occurrence of a singularity seems to be a general feature of the 
simple equation (3.3). It is a consequence of the integrability of the problem 
which is due to the assumption of an isotropic pressure p{r) in (2.6). On 
the other hand the dark contributions =fwo(r) in (2.8) and (2.9) appear with 
different signs. Since the normal matter is probably following the dark in a 
certain way, the pressure components Pj(r) in (2.8) and (2.9) should not be 
the same. We shall investigate this more general situation in the next part 
of this series. 
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